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ABSTRACT 

A  computer  program  assembled  for  the  Hewlett-Packard  9800/ 
Model  20  is  presented  for  calculating  the  rate  of  fire  spread 
according  to  Rothermel's  fire  spread  model. 

OXFORD:  431.6:432:436.  KEYWORDS:  Fire  behavior,  fire 
control,  fire  use,  computer  program,  Hewlett-Packard  9820, 
moisture  of  extinction,  fire  spread  model,  fire  management. 


INTRODUCTION 


I  have  programed  Rothermel ' s •'■  fire  spread  model  for  the  Hewlett  Packard  series 
9800/Model  20  programable  calculator.     This  calculator  quickly  solves  Rothermel 's  set 
of  parametric  equations  for  obtaining  the  rate  of  spread  through  a  porous  fuel  array. 

This  program  was  written  for  two  fuel  categories  with  up  to  nine  size  classes 
within  each  category.     Detailed  data  inputs  given  by  Rothermel-   include  moisture  of 
extinction  (see  Appendix)  for  each  category  and  the  following  fuel  properties  for 
each  size  class  within  a  category: 


•'■Richard  C.  Rothermel,     A  mathematical  model  for  predicting  fire  spread  in  w'  "d- 
land  fuels.     USDA  For.  Serv.  Res.  Pap.  INT-115,  40  p.,  illus.  1972. 


1. 
2. 
3. 
4. 
5. 
6. 
7. 


Load  

Heat  content  (low  heat  value)  

Fractional  moisture  content  

Surface  area-to-volume  ratio  

Total  fractional  mineral  content  

Particle  density  

Silica-free  ash  content  (fractional) 


(Ib./ft.^) 
(B.t .u./lb.) 


(lb. /lb.) 
(ft.-l) 
(lb. /lb.) 
(lb. /ft.  3) 
(lb. /lb.) 
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PROGRAM  OPERATION 


The  fuel  bed  depth  (ft.)  must  be  entered  following  the  fuel  descriptors.  An 
option  is  available  in  the  program  to  weight  the  calculation  of  the  reaction  intensity 
in  each  category  by  the  effective  heating  number,  which  in  turn  is  evaluated  from  the 
characteristic  surface  area-to-volume  ratio  of  the  category.^ 

Finally,  the  slope  and  windspeed  at  midflame  height  (ft./min.)  are  requested  by 
the  programed  calculator.    The  program  is  broken  into  two  parts  because  of  its  large 
size.    After  obtaining  the  rate  of  spread,  the  second  part  of  the  program  can  be 
reinitiated  for  another  rate-of-spread  calculation  with  the  same  fuel  descriptors  and 
a  different  fuel  depth,  windspeed,  and  slope. 

The  program  evaluates  the  following: 

Reaction  intensity  (B.t.u./ft.^-min.) 


Interim  results  can  be  obtained  by  assessing  the  memory  register  (see  Appendix, 
Register  Printout) . 

The  following  program  employs  the  math  ROM  in  block  3;  blocks  1  and  2  are  empty. 
Internal  registers  total  429. 


^William  H.  Frandsen.    Weighting  the  reaction  intensity  by  the  effective  heating 
number.     Intermt.  For.  §  Range  Exp.  Stn.,  USDA  For.  Serv.  Gen.  Tech.  Rep.   (In  prep.). 


Packing  ratio- 
Bulk  density-- 
Rate  of  spread 


(ft.  Vft.3) 
(lb. /ft. 3) 
(ft./min.) 
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Step  Display 


Part  I  Instructions 


ERASE  LOAD  EXECUTE 


3 
4 


0-MANLOD  1-INREG 


5-1      NEXT  MAG  CARDS 


Alternately  insert  magnetic  cards  of  Part  I 
following  the  alphabetical  sequence  and 
EXECUTE  until  NOTE  14  is  no  longer  displayed 
indicating  completion  of  loading. 

If  input  data  are  to  be  entered  by  magnetic 
cards  insert  after  completion  of  the  program 
cards  and  EXECUTE  until  NOTE  14  is  no  longer 
displayed. 

END  RUN  PR0GRA^4 

0^  RUN  PROGRAM  if  data  are  to  be  entered 
manually  through  the  keyboard 

1^  RUN  PROGRAM  if  data  was  entered  from 
magnetic  cards 

Option  1  leads  directly  to  the  end  of  Part  I 
indicated  by  flashing  display  (lOX)  followed 
by  automatic  initiation  of  load  execution. 
Calculator  is  now  ready  for  insertion  of 
Part  II.     Go  to  Part  II. 


5-0      CLASSES  CAT  1=? 


Option  0  begins  manual  entry, 
m  RUN  PROGRAM 

m  =  classes  in  first  category,     m  7^  0 


CLASSES  CAT  2=? 


n  RUN  PROGRAM 

n  =  classes  in  second  category, 
is  only  one  category. 


n  =  0  if  there 


MEXT  1=? 


8  MEXT  2=? 


£  RUN  PROGRAM 

p  =  moisture  of  extinction  for  category  1 
q  RUN  PROGRAM 

q  =  moisture  of  extinction  for  category  2 
ignored  if  n=0  in  step  6. 


LOAD  =  ? 


a.   RUN  PROGRAM 
1  

a^  =  load  (lh./ft.^)  in  class  sequence  i  where  i 

is  sequential  numbering  of  classes  beginning 
with  category  1  through  category  2;  i.e.^  if 
there  are  three  classes  in  both  category  1  and 
category  2  then  i  =  4  for  the  first  class  of 
category  2.     Sequence  number  flashes  between 
displays . 


step 


Display 


Part  I  Instructions 


10  HEAT  CONT=? 


b.   RUN  PROGRAM 
1  

=  heat  content  (B.t.u./lb.)  in  sequence 

class  i. 


11 


FRAC  MOTS  =  ? 


c.   RUN  PROGRAM 
1  

c^  =  fractional  moisture  content  in  sequence 

class  i. 


12 


13 


SA/V  =  ? 


MNRL  CONT  =  ? 


d.   RUN  PROGRAM 
1  

d.  =  surface  area-to-volume  ratio  (ft.  ^)  in 
1 

sequence  class  i. 

e.  RUN  PROGRAM 
1  

e.  =  total  fractional  mineral  content  in 
1 


sequence  class  i 


14 


PART  DENS  =  ? 


f.   RUI^  PROGRAM 
1  — 


f^  =  density  of  fuel  particle  (lb./ft.^)  in 
sequence  class  i. 


15 


SFA  =  ? 


16  NEXT  MAG  CARDS 


q.  RUN  PROGRAM 
^1  ■  

q^  =  silica  free  ash  content  in  sequence 
class  i. 

Steps  9-15  repeated  until  i  =  m  +  n. 

Flashing  display  (lOX)  followed  by  automatic 
initiation  of  load  execution.  Calculator  is 
now  ready  for  insertion  of  Part  II. 
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Step  Display 


Part  II 


Instructions 


1-A 


1-B 


1-C 

2 
3 


FUEL  DEPTH=? 


WT  BY  EPF  HTNG? 


Alternately  insert  program  magnetic  cards  of 
Part  II  following  the  alphabetical  sequence 
and  EXECUTE  until  NOTE  14  is  no  longer  displayed 
indicating  completion  of  loading.     Go  to  Step  2. 

If  Part  II  is  the  beginning  and  Part  I  did  not 
precede  Part  II  in  operation  (see  later 
discussion) ,  then  ERASE  LOAD  EXECUTE . 

Follow  instructions  under  Step  1-A  above  and  go 
to  Step  1-C. 

Insert  magnetic  data  card  and  EXECUTE  until 
NOTE  14  is  no  longer  displayed. 

END  RUN  PROGRAM 

d  RUN  PROGRAM 

d  =  Depth  of  fuel  bed  (ft.). 

Option  to  weight  by  effective  heating  number. 
RUN  PROGRAM 


0-YES  1-NO 


WIND  =  ? 


SLOPE  =  ? 


Option  0,  0_  RUN  PROGRAM 

The  program  will  calculate  the  characteristic 
effective  heating  number,  e. ,  of  category  i 
using  the  characteristic  surface  area-to-volume 
ratio  of  category  i,  replacing  f . ,  the  weight 
factor,  for  category  i  .-^ 

Option  1 ,  j_  RUN  PROGRAM 

Calculation  proceeds  as  described  by  Rothermel 
using  f^,  the  surface  area  weighting  factor, 

w  RUN  PROGRAM 

w  =  windspeed  at  midflame  height  (ft./min.). 

s_  RUN  PROGRAM 

s  =  slope  (Ay/ Ax) . 


If  there  is  only  one  class  in  one  category,  changes  in  fuel  parameters  can  be 
achieved  by  directly  accessing  the  appropriate  memory  register.    However,  if  more  than 
one  class  or  category  is  required  and  changes  are  to  be  made  in  any  of  the  following 
three  parameters:  load,  surface  area-to-volume  ratio,  and  particle  density,  then  it  is 
necessary  to  begin  again  with  Part  I.     One  needs  only  to  enter  changes.     For  all  other 
input  requests,  press  Run  Program. 
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APPENDIX 

Recording  Data 


EXECUTE  THE  FOLLOWING: 

REC  "DA",  R(8(Rl-t-R2)-t-20) 
EXECUTE 


Register  Printout 

THE  FOLLOWING  PROGRAM  LISTS  ALL  REGISTERS  UP  TO 
AND  INCLUDING  REGISTER  R7+16. 


0  -i-  H  ?  7  (  R  1  +  R  2  )  +  3  7 

Bh 

FXIi  U?PRr  H5FXII 
4  •!  F'  R  T   R  R  i!  H  +  1  -i*-  H  ■! 

SPC  Sh 

Ei-ili  h 
1-1 .1  "i  ■•■•i 

>\  -i-  C.  c. 
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Register  Identification 

REGISTERS  QUANTITY 

0  0  -  Manual  load,  1  -  data  entered  by  magnetic  card 

1  Number  of  classes  in  category  1 

2  Number  of  classes  in  category  2 

3  Characteristic  SA/V  of  category  1  (SA/V  =  Particle 

surface  area-to-volume  .ratio) 

4  Characteristic  SA/V  of  category  2 

5  Characteristic  SA/V  of  entire  fuel  bed 

6  Unused 

7  7(R1+R2)+21 

8  f  ^  ,  weighting  parameter  for  category  1 

9  f^,  weighting  parameter  for  category  2 

10  Fuel  depth 

11  Corrected  organic  load  of  category  1 

12  Corrected  organic  load  of  category  2 

13  Weighted  heat  content  of  category  1 

14  Weighted  heat  content  of  category  2 

15  Mineral  damping  coefficient  of  category  1 

16  Mineral  damping  coefficient  of  category  2 

17  Moisture  damping  coefficient  of  category  1 

18  Moisture  damping  coefficient  of  category  2 

19  Moisture  of  extinction  of  category  1 

20  Moisture  of  extinction  of  category  2 
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REGISTERS 
21,28. . 
22,29. . 
23,30. . 
24,31. . 
25,32. . 
26,33. . 
27,34. . 
+ 

R7-1 
R7 


8(Rl+R2)+20 

R7  +  10 

R7+11 

R7+12 

R7+13 

R7+14 

R7  +  15 

R7+16 


QUANTITY 

Load  (i  J  j  ) 

Heat  content  (i^j) 

Fractional  moisture  content  (ijj) 

SA/V  Ci,j) 

Total  fractional  mineral  content  (ijj) 

Particle  density  (ijj) 

Silica-free  ash  content  (ijj) 

j  =  1  (category  1)  i  goes  from  1  to  Rl 
j  =  2  (category  2)  i  goes  from  1  to  R2 

=  7(Rl+R23+20  last  of  input  data  list 

Weighting  parameters  f(l,l) 

4- 

f(l,Rl) 

f(2,l) 

i 

f(2,R2) 

R/B  6      =  optimum  packing  ratio 

op  op         r  r  a 

Wind 

Slope 

I  =  reaction  intensity  (B. t .u./ft . ^-min.) 

K 

=  propagating  flux  (B . t . u. /ft . ^-min . ) 
B  =  packing  ratio 

=  bulk  density  (lb. /ft. 3) 
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Program  Listing 


FIRE  SPBEAD  MODEL 
Part  I 


1  S ' 

7.';R1+R2i  -21 -RTF- 
lb : 

R  (  H  7  +  2  1)  R  I  V?  +  4  ! 
/R  i:R7  +  26  i  -:^R  I  R?  +  H 
::  R 


FIRE  SPREAD  MODEL 
Part  II 


Efir    0  riRHLon  i- 

1 KREG" ^  RQi  IF  R0  = 
IjGTO  "STRT"!- 
1 : 

ENT  "CL.R3SES  CRT 
i='?"  1  Rl  J  "CLASSES 

C:RT2  =  '?  • )  R2h 
2 ! 

0-  *R;EHT   "NEXT   1  = 
•?",R19;  IF  R2>o;- 
ENT   "NEXT  2  =  ?"!-R 

2ei;GSB  "r'h 

■~l  ■ 

EHT   "LORri  =  '?"  J  R  (R 
7+21) ;GSB  "l"h 
4: 

ENT   "HERT  C:ONT  =  '> 
"  J  R  f.fl7  +  22:' ;  GSB  " 
1  "h 
5: 

EHT   "FRflC  Mi:iIS  =  ':^ 
"  J  R  i:R7  +  23  .i  j  GSB  " 
1  '> 
6 : 

ENT  "SR/V=?" >R (R 
7  +  24);GSB  "r'h 

ENT  •riNRL  i::ONT  =  ';' 
"  J  R  i:R7  +  35  3  ;  GSB  " 
1  "h 

1-  i ' 

ENT   "PART  IiENS  =  ":' 
" J  R  i  R7  +  2b )  ; GSB  " 
1  "h 
9 ' 

ENT   "SFA  =  '?"  ?  R  i  A7 

+  2  7  j  h 

10: 

GTO  +21- 
1  1  : 

" 1" ; r+i^b; Fxn 

nSP  BjDSP  BMiSP 
BjDSP  BjDSP  E; 
RET  h 
12? 

IF  ri+R2=b;gto  + 

2h 
13: 

b-*r;gto  31- 

14: 

"strt" ; 0^Ah 


RHl^R j IF  Rl+R2>fi 

;gti:i  --if- 

18: 

0  ->  R  4  B  ->  C I-- 
.19:  ■ 

IF  R1>R5  B  +  R  i:R7  +  H 

!! -^b; r+i^r; JMP  ei- 

29: 

IF  R2>0; IF  R1+R2 

>  R ;  c  -r  R  ( R  7 + R  :i  ^  c ;  R 

+  l^AjjriP  HiF 
21  : 

e->RF 

IF  R  i  >  R  ?  R  i:  R  7  +  fl  j  .-■ 
B^R  i;  R7  +  R  j  ?  R+  1  -^R  j 
JMP  ©F 

IF  R2>0; IF  R1+R2 
>R:  R  rR7  +  R:!  /C^R  (R 

7+r:!  ;  r  +  i^r;  jmp  e 
I- 

24: 

B-  (B  +  C) ^RS;  IF  R2 

>  0  ?  C (  B  +  C  j  -^R'ih 

III  ••'  " 

0  -i^  fl   B  -i^  C  F 
26 : 

IF  R1>R? B+R (R7+R 
j  R  CR7  +  21 j / ( 1+R (fi 
7  +  25  3  j  -^-BF 

IF  R"l  >  R  ?  R+ 1  ii 
GTO  -IF 
2'  o : 

IF  R2>05 IF  R1+R2 
>R5  C  +  R  (R7  +  R:i  R  i:R7 
+2ij/!l+RiR7+25) 

j  HHJl- 

29 : 

IF  R2>0? IF  R1+R2 
>fl;R+l+H;GTO  -IF 
30 : 

B-^R  1  1  ;  C-^R  1  2F 

:>i  i 

0^RF 

DSP   "HE>rr  MRG  CR 

Rjjs" ; r  +  i-i^r;  IF  f\i 
i0;jriF  0F 

END    SLGD  F 
R30S 


0 : 

0^ri-^B-K:F 
1  : 

IF   R1>R? B+R f R"+H 

.1 R  I  R7+22:i  ->b;  R+1-^ 
r;jnp  0F 

IF  R2>0.5  IF  R1+R2 
>R;  C  +  R  (R7  +  Ri  R  i:R7 

+22  j  -^c ;  H+ 1  -^A ; 
jrip  0F 

o : 

B-^RlS;  C^R14F 
4 : 

0->R->B-^C:h 

■J  « 

iF  ri>r; B+R i;R7+R 
)  R  I  R  7  +  2  7  i  -i^  B ;  R  +  1  -> 
RjJMP  0F 
6 : 

IF  R2>0; IF  R1+R2 
>  R  ;  C  +  R  I  R  7  +  R  3  R  (  fl  7 
+  27  )  H^C  ;  R+  1  ->R  j 
JMP  0F 

.  174Bt  (-.  19)  -^RIS 
IIF  R2>0;  .  174i::t  i 
-,193  ^RlbF 

0  ■ 

g-i^R-i^D^CF 

9 : 

IF  ri>r;  B+R  i:R7h  H 

1  R  C  R  7  +  2  3  3  ^  B  °  H  +  1 
R5.JMP  0F 

10: 

IF  R2>0; IF  R1+R2 
>R; C  +  R (R7  +  R3  R (R" 
+  23  3  ^D;  R+I-^r: 
...IMP  0F 

1   1  5 

e/R19^Z»  GSB  MIi^ 
F 

12: 

Z  +  R  1  7  F 
1  3 ' 

IF  R2>0; C/R20-Z? 

GSB  'MIi'-F 

14: 

Z  +  R  1  S  5  G  I  0  +  3  F 
15: 

■^MIi"  ?  IF  Zi  1  :  1-2. 
59Z  +  5»  1  lZt2"3f 52 
ZtS^ZjRET  F 
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16: 

0^Z;R£T  h 

175 

18: 

IF  R1>R!  B  +  R  i:R7  +  R 

n;  jMP  01- 

19: 

IF  R2>0; IF  R1+R2 
>  H  ?  C  +  R  i:  R  7 1-  fl  )  R  i:  R 

2  4  )  ^  C  j  fl  +  1  ^  R  ; 
..IMP  01- 
28 : 

B  ^  R  3 ;  i:  ^  F;  4 ;  r  s  r  3  +  f: 

9R4^R5I- 

21  s 

ENT    TUEL   riEPTH  = 
■?  "  ?  R  1  U  5  FKIi   3  j 
SPC   2  5PRT  "FUEL 
riEPTH  =  "  5  R  1  0  ;  0->R-^ 
Bh 

IF  R1+R2>R; B  +  R 'R 
7  +  2  1  j    R  C  R  7  +  2  6  j  E; 

;R+i^R;.jNP  01- 

B/R  1  ©T^B  5  3  .  34SR5t 

i -.31 89  j  H^:>:: ;  B/x-?x 

;  B^R  i:R7+15)I- 
24: 

i/ (4, 774R5t. 1-7. 
27)  -^R;  R5-t  l  .  5/  (49 
5  +  .  0  5  9  4  R  5 1 1  .  5  j  -i*  C 

h 

25 : 

C>::tREXP   (R  i  i-x:i  i 

■^c;  x^R  i;r7  +  10:i  h 

2  6 : 

ENT  "WT  BY  EFF  H 
THG'^"  J  n?  "0-YES  1 
-HO"  J  Rh 

1^   i'  ■ 

IF  R=l  ;  RSh^X!  R9^Y 

;gtij  +21- 

E;<P  (-138/R3) 
IF  R2>0;EXP  (-13 
8/R4.1  ^Yh 
29 : 

CXRl  lR13R15Rl'7-^X 
5  CYR12R14RlbR18-* 

YF 
30 : 

SPC  ZjFxd  6;prt 

"RCTN  INT-l  =  "j>:; 
IF  R2>0;PRT  "RCT 
N  INT~2=^"?Yh 


31  : 

X  +  Y-i^Y-^R  (R7+13;i  ; 
PRT   "TOT  RCTN  IN 
T=" ?  Yh 

YEXP  i:  (  .  792+.  631 
R5t.  5  j  i;b+.  1  :i  )  /  I  1 
9  2  +  .  2  5  9  5  R  5  )  R  (  R  ?' 
+  1  4  j  h 

ENT   "I..IINri  =  '"'"  !.  R  (R 
7+i  1  j  1.  "SLOPE  =  '"'"  J 
R  i  R7+12:i  ;  SPC  21- 
34  : 

PPT    '■MIHD="  1  R  i:R7 
+  1  1  j  r  "SLOPE=" J  R ( 
R?+12:i  h 
35 : 

7. 47EXP    (-. 133R5 
•t-  „  5  5  )  ^  C  ;  .0252  6  R  5 
t.  54->Bh 
36 : 

.  715EXP  i:-3.  59E- 
4R5:i  -^Rh 

•~i  ^  ■ 

._t  1  a 

5 .  '2  7  5  R  (  R  7  +  1  5  i  t  (  - 
.  3  j  R  i  R7+  1  2  j  t2-*YF 

CR  i:  R7+  1  1  .1  tBR  (  R7  + 
1  0  i  t  (  -R  j  H^XF 
3  9 : 

R  I  R  7  +  1  4  )  i:  1  +  X  +  Y  j  ^ 

H  ;  0  R  't'  f 

40: 

IF    .  9R  i:R7+13:i  iP  ( 
R7+1  1  i  5  SPC    I  ; 
PRT    ".9RCTN  INTi 
UINU"  ;  0-^fi^YF 
41  : 

IF   P1+R2>R; Y+R (7 
R  +  2  1  )    Y  !!  R  +  1   R  5 
-IMP  0h 
42: 

Y  /  R  1  0  ->  Y  ^  R  (  R  7  +  1  6  ) 
;  0   R  -i^  B  ^  C  F 

IF  RDflJGSB  "QG" 
F 

44: 

IF  R1>R5 B+R (R7+R 
j  Z-^Bj  R+l-^Rj  GTO  - 
IF 
45: 

IF  R2>0; IF  R1+R2 
>R;GSB  "QG'> 


46 : 

IF  R2>05 IF  R1+R2 
>R;C  +  Ri:R7  +  R)2->Cj 
R+l-fflJGTO  -Ih 
47: 

Y  i:br8+cr9)  ^y; 

GTO  +2h 
48 : 

"QG";EXP  (-138/R 
i:R7  +  24:i  )  (250+1  1  1 
6R  (R7  +  23)  j  -^ZS 
RET  F 
49: 

FXn  6; SPC  2; PRT 
'PRCKING  RRTIO=" 
1  R  (R7+15  j  ■>  "BULK 
riENSITY="  J  R  (R7+1 

6:>  F 
50: 

SPC  2;FXD  3; PRT 

"RRTE  Of  SPRERD= 

"  J  X/YF 

51  : 

END  F 

R203 
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Sample  Register  Printout 


fUMfAC  it>MO     y  .  y  u  0  U 

1 

eOiSfS  lAj  CAT  I  4.00  0  0 
CU»fi  /A>  Ca-t  2-   4  .  0  0  0  0 


(T  CAT  / 

cr  corz. 


i  2  4  S  .  0 0  3 


-J  1.^  ■. 


9  4 7 .249  6 

0  .  0  0  0  0 


i 


I*' 


0  0  0  0 

5S4  6 

Q 

4  i  54 

10 
0  0  0  0 
1  1 
015? 
1  2 

0  3  6  3 
13 

k,    8955.  91  38 

14 

8  984.40  3  2 
(7J,  .3501 
(l*)^  .3b 6  4 
.6503 
(T-),.7886 

(Mtot)  3.1300 

20 
.2000 


.00. 


"1  r,  f;,  .  000  \  \ 


7  5  U  0 

V  J, 


'  ©"fcL  2  000  .  0tiOti 

(&;K;  -  0  550 

26 

CPp),j    20.  0  000 

CSe)ii  .0350 


7^/ 


025fc 

2  9 

0  0  0  0 

3  0 
5500 

31 
0  y  0  0 

0150 

0  0  0  0 
34 

0  150 


y  y  y  y 

■J"  1' 

5500 


(0^,  .010 

^\  3327 

(4  133 
3  1 


y  0  y  y 

3  9 
015  0 

4  0 

0  0  0  0 
4  1 

0150 


4  2 
.  0043 
4  3 

.  0000 

4  4 
.550  0 

45 

6  7 , 

.  0000 

4  6 
.0150 

47 

4  7 

.  0000 

4  8 
.0150 

4 

49 

.  0002 

50 

895  3 

.      0  0  0 

51 

.  0200 
52 

r  -I 
•J 

.  0150 

54 

2  0  .  0  0  0  0 
.0150 


56 

57 

58 

.  02  00 

5  9 

500. 0  0  0  0 
60 

.  0150 

61 

1  8  .  0  0  0  0 

.0150 

6  3 
.  00  33 
6  4 

8  32  7. 0  0  0  0 

65 

.  0  200 

66 

1  5  3  .  0  0  0  0 

67 

.0150 

6  S 

4  6.0  0  0  0 

69 

.0150 

/-I 


7  1 

3  9  3  .  0  0  0  0 
.  0200 
73. 0000 

(  t 
.015  0 

75 

51.000  0 
76 
.  0  1  50 


77 

.5113 

78 

I. 

.  4540 

79 

.  0307 
80 

.  0040 

81 

.0184 

82 

.  971  0 

33 

.0101 

84 

.  0004 
85 

0 

.  0000 

0 

86 

.  0000 

87 

U)JA<D      17  6.  0  0  0  0 

8  9 

SLOPE.  .4100 

90 

9  1  9  ,  0  2  0  0 
■91 

X,  19.3955 
92 

^  .0022 
fs  .0452 
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Sample  Output 


A  sample  output  using  the  preceding  sample  register  data  is  sho;vn  below. 

FUEL  IiEPTH  = 

2 .  688 

RCTH  INT-i= 

256. 233431 

TOT   RCTN  IHT= 

9  1  9  .  6  2  6  6  6  ? 

WINri  = 

1  r  6  .  0  0  6  y  H  fi 

SLuPE= 

.410  e  0  0 

p  H  c  k:  i\\  g  k  r  t  1 0 = 

=002224 
BULK  IiEHSITV  = 

.045150 


RhTE  of  SPREhI:  = 
13. 292 


Moisture  of  Extinction 

The  following  program  (written  for  the  HP  9820)  will  evaluate  the  moisture  of 
extinction  necessary  as  an  input  to  the  fire  spread  model. 

Three  inputs  are  requested:   (1)  dead  fine  fuel,   (2)  live  fine  fuel,  and  (3)  dead 
moisture  content.    The  first  two  are  fuel  loads  and  are  entered  in  Ib./ft.^.  The 
operator  can  decide  which  of  the  fine  fuels  are  contributing  to  the  desiccation  of  the 
live  fine  fuels  and  which  of  the  live  fuels  are  affected  by  the  desiccation  process. 
As  a  rule,  the  size  classes  considered  to  be  fine  fuels  and  used  to  calculate  the  mois- 
ture of  extinction  are:  foliage,  0-inch  to  %-inch,  and^-inch  to    -inch  from  the  dead 
fuel;  and  foliage  and  0-inch  to  %-inch  from  the  live  fuel.     (Note  that  foliage  is  a 
separate  class.)     The  dead  moisture  content  is  taken  from  the  0-inch  to  ?3-inch  dead 
fuel  size  class.     (Moisture  contents  of  the  dead  fine  fuel  size  classes  should  not 
differ  greatly.)     Outputs  include  the  ratio  of  the  live  fine  fuel  load  to  the  sum  of 
both  the  live  and  dead  fine  fuel  loads  and  the  moisture  of  extinction. 
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ENT  "liERD  FINE  F 
UEL';:"'  H  fl?  "LIVE  FI 
HE  FUEL?"hE:|- 

ENT  ■■HERD  NO  IS  C 
uNT""'"  C- pRT  "DEFi 
D  FINE  FIJEL="?R; 
B .  -  •  i  R  +  B  :i  -i^  R  F 


PRT  "LIVE  FINE  F 
UEL="  11  "DERD  NO 
I  S   CONT-  '■  ?  C  j  SRC 


PRT  ■■LIVE  TO  TOT 
RL  '  "FINE  FUEL  R 
h  J  i  iJ    '■  ?  R  F 

2  ::  9  (  1  Fl    {  1  -  1.  0  C  '  3 

:  .  ■'  !■■!  „  2  2  6:  Z. 

SPC    i;pRT    "EKT IN 
CTIGN"  .  "MOISTURE 
CGNT=^'  5  ZF 


ENIi  F 
R  3  9  8 


Sample  Output 

DERIi   FINE   FUEL  = 

.  0 

LIVE   FINE  FUEL= 

.  0 

D E I ■•! l!   N !J  I  S   C \J\\J  = 


LIVE  Tn__TOTRL 
r  J.  i'  -I  i:::    i  "  u  b.  L    R  R  i  i  L 


ErtTINCTIGN 

Ml!  I  GTURE  CONT- 

d 
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